#ifndef __MKL_UTILS_NURBS__
#define __MKL_UTILS_NURBS__

#ifdef __MKL__
namespace MKL_wrapper {

void SolvePardiso(const boost::numeric::ublas::mapped_matrix<double>& A, const boost::numeric::ublas::vector<double, std::vector<double> >& b, boost::numeric::ublas::vector<double, std::vector<double> >& x, const bool symmetric);
void SolvePardiso(const boost::numeric::ublas::mapped_matrix<double>& A, const std::vector< boost::numeric::ublas::vector<double, std::vector<double> > >& b, std::vector< boost::numeric::ublas::vector<double, std::vector<double> > >& x, const bool symmetric);

}
#endif //__MKL__

#endif //__MKL_UTILS_NURBS__
